Theory of molecular crowding in Brownian hard-sphere liquids 
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We derive an analytical pair potential of mean force for Brownian molecules in the liquid-state. 
Our approach accounts for many-particle correlations of crowding particles of the liquid, and for 
diffusive transport across the spatially modulated local density of crowders in the dense environ- 
ment. Specializing on the limit of equal-size particles, we show that this diffusive transport leads 
to additional density- and structure-dependent terms in the interaction potential, and to a much 
stronger attraction (by a factor ~ 4 at average volume fraction of crowders (j>o = 0.25) than in the 
standard depletion interaction where the diffusive effects are neglected. As an illustration of the 
theory, we use it to study the size of a polymer chain in a solution of inert crowders. Even in the 
case of athermal background solvent, when a classical chain should be fully swollen, we find a sharp 
coil-globule transition of the ideal chain collapsing at a critical value of the crowder volume fraction 
4> c « 0.145. 

PACS numbers: 82.20.Uv, 34.10. +x, 95.30.Qd 



I. INTRODUCTION 

The theory of the liquid state relies heavily on the sta- 
tistical mechanics of hard-sphere systems [l| . One of the 
most remarkable among these results is the Percus-Yevick 
integral equation, the solution of which yields an accu- 
rate description of the structure of hard-sphere liquids 
in terms of the radial distribution function g(r). The 
latter is related to the probability of finding a second 
particle at a position r along the outward radial coordi- 
nate measured from the test particle. g(r) typically fea- 
tures a peak near contact (the first coordination shell) 
associated with an enhanced probability of finding a sec- 
ond particle near contact. This is due to the osmotic 
pressure exerted by all the other particles in the system 
which is not balanced in the gap between the two parti- 
cles near contact, an effect which is known as depletion 
attraction. This effect is clearly the result of collective 
dynamics ( "more is different" ) since it introduces an at- 
tractive interaction which depends on the density, and is 
not present originally in the pair interaction between two 
hard particles. In ensembles of hard spheres of equal size, 
this attraction induced by many-body effects is known to 
be weak: even at high density it is barely of the order of 
the thermal energy scale kT Q. This attraction, how- 
ever, may become substantial in binary mixtures of hard 
spheres where spheres of two different sizes coexist: this 
is an important attractive interaction that determines 
the phase behavior of colloid-polymer mixtures [H and 
its first theoretical description (valid for dilute systems) 
was given long ago in a famous work by Asakura and 
Oosawa Q. This attraction plays an important role in 
biological systems where, e.g., it is a factor controlling 
cellular organization 

However, in both colloidal and biological systems, the 
particles are not simply hard spheres subject to New- 
tonian dynamics. Instead, they are Brownian particles 
obeying Langevin dynamics and they are sensitive to the 
local chemical potential of solvent [fj. Since this aspect 



has been neglected both in theoretical and simulation 
studies of depletion in the past, it is an interesting ques- 
tion to ask: what is the effect of the dense solvent and 
the diffusive nature of particle dynamics on the depletion 
attraction. Here we show that many-particle correlations 
that lead to the modulated radial distribution g(r) can 
strongly affect the local osmotic pressure profile between 
two test particles. As a result, diffusive driving forces 
arise which push the two particles together more strongly 
than if the dynamics were purely Newtonian, as has been 
assumed in [4j and ever since. This is reflected in a much 
stronger depletion attraction for Brownian particles as 
compared to hard spheres in vacuo, such that the at- 
traction minimum can be substantial also in systems of 
equal-size particles. 

We then apply this theory to calculate the coil-globule 
transition of ideal polymer chains in solution with Brow- 
nian particles of approximately the same size as the 
monomers. It was observed 40 years ago that adding 
polyethylene glycol (PEG) molecules to solutions of DNA 
in water leads to the collapse or condensation of DNA as 
soon as a critical concentration of PEG is reached 0. 
Several theories have been proposed to explain this 
"crowding effect" , which has widespread implications as 
the biological function of all biopolymers is strongly de- 
pendent on their size and conformation, in particular, 
on whether they are in a swollen or collapsed state 
This stark distinction in biological function is reflected 
in the names of their "denatured" and "native" states of 
biopolymer molecules. One of the most familiar cases is 
the packing of DNA, in dense cellular environment, and 
in vitro experiments where the crowding was artificially 
tested [§, [l0] • In proteins subtle changes of molecular en- 
vironment, instead of natural folding, could trigger their 
aggregation into toxic assemblies (such as the amyloid 
fibrils responsible for neuro-degenerative diseases) pd| . 
It has been shown that folding can be promoted by the 
addition of inert crowders via a mechanism which bears 
many similarities with the crowding-induced globular col- 
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lapse of generic polymers fl2| . 

The effective interaction between pairs of polymer seg- 
ments in the presence of quenched randomly distributed 
impurities is a classical problem in polymer physics, 
originally solved by Edwards, Muthukumar and Baum- 
gaertner [IH [l4[ . When the impurities are mobile (an- 
nealed), a similar attraction also exists and is likewise 
due to entropic effects, this time in the form of depletion 
forces [Hj]. Despite the similarities, the mobile crowding 
problem has a different formulation because the effective 
potential in this case has to be related to the whole diffu- 
sive processes acting between the polymer segments and 
the Brownian crowders. The interplay between diffusion 
in a dense environment with a liquid structure and the 
depletion effect, is the focus of our work. 

The theories proposed to explain this phenomenon at- 
tempt to account for the many-particle correlations due 
to crowders, which induce an effective attraction between 
pairs of monomers in the chain [l6l - |2]| . This is reflected 
in a decreasing second virial coefficient upon increasing 
the volume fraction of crowders in the system, until the 
second virial coefficient finally changes sign, from posi- 
tive to negative, and a transition from coil to globule of 
the polymer occurs [2(|. None of these theories, how- 
ever, account for the role of the microscopic liquid struc- 
ture in this process and they all generically attribute the 
increased effective attraction to the increased global (uni- 
form) osmotic pressure exerted by the crowding particles. 
These approaches also neglect the role of the diffusive 
dynamics. In contrast, our approach uses as its input 
the radial distribution function of the fluid <?(r), which 
develops a modulated structure at increased densities. 
It thus allows us to obtain additional forces generated 
by diffusive transport across local gradients of the os- 
motic pressure (and the corresponding chemical poten- 
tial). At higher densities this effect turns out to be much 
stronger than the standard depletion forces, and it of 
course vanishes in the low-density limit where one recov- 
ers the standard depletion attraction given, for equal size 
spheres, by — \ng(r). The fact that all depletion-based 
theories underestimate the crowding effect manifests it- 
self, for instance, in that the critical density of crowders 
for the polymer collapse, derived in (l9j . diverges to in- 
finity when the size of the crowder particle matches the 
size of the polymer Kuhn length - clearly an unphysical 
outcome. 

There are several steps in our development: (1) We 
first derive the analytical form of this effective pair poten- 
tial between monomers, by accounting for multi-particle 
correlations and also for the presence of the solvent, 
which shows an attractive well, increasing at higher con- 
centration of crowders. (2) This effective potential re- 
quires an expression for the radial distribution function 
g(r), which we calculate using the currently most accu- 
rate analytical approach based on the Wertheim solution 
of the Percus-Yevick equations. (3) Using this effective 
potential we then calculate the second virial coefficient 
for monomer- monomer interaction of a polymer chain in 



good solvent, and use a simple scheme of Flory mean-field 
free energy to qualitatively describe the chain collapse at 
a critical volume fraction of crowders. Our theory uses 
only one serious assumption, the local density approxi- 
mation (LDA), the validity of which is carefully exam- 
ined below. Therefore, this analysis can be applied to 
any kind of system, provided that the structure factor of 
the monomer-crowder mixture can be extracted from e.g. 
radiation scattering experiments. 



II. DERIVATION 

Consider a Brownian particle moving in a sea of other 
similar Brownian particles, which all mutually interact 
via excluded volume (purely hard-sphere particles). We 
let the particles be embedded in a background solvent, 
which is a simple liquid providing the thermal bath. Let 
us take some arbitrarily chosen point of the system as the 
origin and be r the distance vector of the tagged parti- 
cle from the origin. In the local density approximation 
(LDA) [l|, let Hs(r) and ^e(r) be the chemical potentials 
of the solvent and of the Brownian particle, respectively. 
Then the local driving forces acting on a solvent molecule 
and on a Brownian particle are given by F s = — V/z s (r) 
and Fb = — V/is(r), respectively. By its definition, 
the spatial variation of chemical potential is related to 
the variation in the free density per particle associated 
with the fact that both the Brownian particles and the 
background solvent molecules have a locally higher free 
energy in regions where the concentration of Brownian 
molecules or, respectively, solvent molecules, is higher. 
As a consequence a driving force caused by the gradient 
of chemical potential acts to bring Brownian molecules 
from region where their concentration is higher into re- 
gions where their concentration is lower, and a similar 
driving force acts on the solvent molecules. At steady 
state, the two forces must be separately equal to zero. 
In fact, the two chemical potentials at constant pressure 
and temperature are related by the Gibbs-Duhcm equa- 
tion: p{r)V ^s(r) + p s (r)\7 fj, s (r) = 0, where p(r) is the 
local density of Brownian particles and Psir) is the local 
density of the solvent. Further, the mass conservation re- 
lation is valid: VBp(r)+v s p s (r) = I, where vb and v s are 
the volumes of the Brownian particle and of the solvent 
molecule, respectively. The relevant, net force which acts 
on the Brownian particle is the one relative to the solvent, 
which is obtained by subtracting the force locally experi- 
enced by the background solvent. In order to ensure that 
the unit volume of background solvent is at rest with re- 
spect to the particle, we need to consider forces per unit 
volume instead of the bare forces @. Hence we have: 
F/vb = Fb/vb — F s /v s , and then the net force acting 
on a single particle is given by: F = Fb — (vb/v s )F s . 
Substituting v s and p s , after a simple algebra we obtain 
the following general expression for the net force acting 
on a Brownian particle moving in the solvent populated 
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with other Brownian particles: 
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■Vn(p(r)), 



(1) 



where the local osmotic pressure is defined as: H(p(r)) = 
— (p s (r)—p,®)/v s , with p® the constant chemical potential 
of the pure solvent that does not contain any Brownian 
particles. This expression for the force is a consequence 
of our using the LDA and is only valid when the local 
variation in density is smooth, that is Vp/p <C l/£ with 
£ the characteristic range of oscillations of the radial dis- 
tribution function g(r) plU^. It is clear, e.g. from the 
Fig. [3 that the interesting effects of crowding-induced 
attraction may occur at sufficiently low densities where 
this LDA condition is satisfied. 

Equation (TTJ) can also be derived from dimensional ar- 
guments. In a suspension of Brownian particles, n is the 
pressure exerted by the suspension itself. Then —VII 
represents the local density of force acting in the sus- 
pension [iH . If p is the number density of particles, it 
follows that the force acting on average upon a single par- 
ticle of the suspension is given by: — (l/p)VIL This force 
is zero whenever the system is spatially homogeneous, as 
it ought to be, e.g. for infinitely diluted suspensions in 
the ideal-gas regime. If, however, the distribution of the 
particles is not homogeneous, this diffusive force acts to 
spread out local concentration gradients. As a result of 
this force, the particle acquires a velocity which points 
along the direction of local steepest decrease of the os- 
motic pressure. The close connection with diffusive mo- 
tion becomes evident upon associating it with the diffu- 
sive flux j = pu, where the overdamped particle velocity 
is: u = (D/kT)F. By applying the continuity opera- 
tor to this flux: p = — V • j, in the low-density limit of 
II = kT p, one recovers the Fickian diffusion equation: 
p = DV 2 p. Within LDA, one replaces 



vn = Vp 



an 



(2) 



thus obtaining the diffusion equation valid for non-dilute 
and locally inhomogencous systems: 



dp(r) 
dt 



D c (p)V 2 p(r). 



(3) 



The collective diffusion coefficient D c = (D/kT)dH/dp 
takes care of the increased density effects [6| ■ This deriva- 
tion highlights the fact that Eq.([T]) represents the driving 
force acting on a Brownian particle in non-dilute and lo- 
cally inhomogeneous systems, and is consistent with the 
particle migration being governed by diffusion, which in 
turn is driven by gradients of the chemical potential. 

Let us now consider that the center of frame from 
which r is measured coincides with a test particle and 
consider a second particle approaching the test particle, 
in the moving spherical coordinates frame centered on 
the test particle. The work required to bring the second 



particle to a point r measured in the moving frame of 
the first particle is given by 



W = — Fdr = 



1 



p{r) dr 



U{p(r))dr. (4) 



In the second equality we used Eq.(fT]) and the fact that 
with hard-spheres the distribution of particles is spa- 
tially isotropic, thus only the radial distance r = \r\ 
between the two monomers is a variable of the effective 
pair-potential energy V e g(r) = W acting between the 
two Brownian monomers. To evaluate V e g (r) we need an 
equation of state for the crowders to relate the osmotic 
pressure and the density, and we will assume that this 
relation also applies locally which, once more, amounts 
to making use of LDA. For the equation of state we can 
use the Carnahan-Starling (CS) equation of state which 
is accurate up to quite high liquid densities [l[ and we 
write it as: 



1 n(p(r)) I + (f)(r) + (f>(r) 2 — (j>(r) 3 



kT p(r) 



(l-#r))* 



(5) 



where 4>{r) = (ira 3 /6) • p(r) is the volume (packing) frac- 
tion of crowders (a chosen as the diameter of the Brow- 
nian particle). The CS equation of state is frequently 
used in standard approaches for the calculation of de- 
pletion forces. However, we should remark that in those 
cases the aim is to account for density effects on the free 
energy of interaction which for Newtonian hard-spheres 
is a purely mechanical effect. That is, one takes the work 
W oc n • Va, with Vd is the depletion volume (a function 
of r), and the CS equation is used to relate the global 
(r-independent) osmotic pressure of the solution n to the 
overall density of the crowders <po- In those approaches 
II and 4> are not allowed to vary with r and therefore the 
driving force in Eq.(fT]) due to local variation in particle 
density is neglected. See Fig. Q] for a visual explanation 
of this contrast. 

Using this CS expression ([5]) in the r.h.s. of Eq. (fj| 
and integrating we obtain 



Vcsjr) 
kT 



fo(r) -1)3 (</> -l) 3 



In- 



(r) 



(6) 



where it is important to distinguish the local packing 
fraction of crowders (j>(r) from the mean fraction 0o a t 
r — > oo. By construction, this is an effective pair poten- 
tial of interaction between two Brownian particles. In 
addition to the unbalanced osmotic pressure in the de- 
pletion gap near contact (which is expressed by the stan- 
dard — \a\<p(r) / 4>vi\ = — hig(r) term), it accounts also for 
the diffusive drift pushing the two particles against each 
other because of the local liquid structure, expressed by 
the two fractions in Eq.®. 

Equation ^ is the key result of this work. Note that 
the — lng(r) term represents the depletion effect, which 
many of the earlier important publications have inves- 
tigated thoroughly 0, [H, [13 ■ For dilute systems wc 
have that the radial density profile around the monomer 
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FIG. 1: (Color online) A: The classical depletion effect when 
no third particle can be accommodated in the gap between 
the two tagged particles. As a consequence, the particles are 
pushed towards each other by a force proportional to the 
global osmotic pressure of the suspension II which remains 
unbalanced in the gap. B: The generic effect of the diffusive 
driving force, Eq.fTJ, proportional to the local gradient of the 
osmotic pressure and driving the particle down the concen- 
tration gradient. This force must also be active on the two 
test particles of panel A, because upon moving closer together 
they are effectively moving across a local concentration gra- 
dient of a dense fluid of crowders. 
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FIG. 2: (Color online) Radial distribution function and effec- 
tive pair potential for <f>o = 0.25 calculated using the Percus- 
Yevick theory [2|| . There are no free parameters in this ather- 
mal system with only excluded volume between particles in 
contact. The dashed line shows the prediction of the depletion 
theory, — lng(r), for comparison. 



tends to become flat, </>(r) — > 4>q for any r. In this 
limit, one can readily check via de L'Hopital's rule that: 
lim^o V e s = —kTlag(r). This limit is a well- known 
law of statistical mechanics [l[ which holds in the di- 
lute ideal-gas limit independently of the particular form 
of the pair-potential. From another perspective, since 
this theory describes the effective pair potential between 
particles of a "Brownian fluid" , this result means that in 
the dilute limit it does not matter whether the particles 
are Brownian or not, which is an important observation 
on its own. Further, our result extends the fundamental 
potential of mean force from the standard simple fluids 
to the Brownian fluids. 

As we have explained above, depletion alone under- 
estimates the effective attraction in Brownian systems, 
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FIG. 3: (Color online) The effective potential upon increas- 
ing the mean volume fraction of Brownian crowders, cj>o = 
0.1, 0.15, 0.2 and 0.25 for the successive curves. 



which is explicitly illustrated by plotting the two forms 
of effective potential in Fig. [3J after substituting the ex- 
plicit form of g(r) at a given mean density 4>o as explained 
in the Appendix [26| . This effective pair potential is plot- 
ted for cf>o = 0.25, which is a relatively high density still 
well within the fluid region of the phase diagram of hard 
spheres. It is evident that the effective potential devel- 
ops a short-range attractive well aligned with the first 
coordination shell of the radial distribution function. It 
is also clear that the attraction minimum is significantly 
deeper, by a factor ~ 4, than the one given by the stan- 
dard depletion potential — In g(r) alone which suffices to 
describe hard spheres in vacuo (that is, ignoring the dif- 
fusion of particles in the suspending matrix). The two 
new terms in the r.h.s. of Eq.© turn out to be a major 
effect not captured by earlier theories. 

Figure [3] plots the evolution of the effective potential 
V e f[ (r) on increasing the mean concentration of crowders 
in solution, illustrating that the depth of the attractive 
well exceeds kT at (j>o > 0.15. It is interesting to note 
that also the local maximum, the energy barrier in cor- 
respondence of the local minimum in the g(r), is signifi- 
cantly larger for Brownian particles than for Newtonian 
hard spheres. This can be understood by recalling that 
the diffusive driving force acts to push the particles down 
the p and II gradients. Upon approaching each other the 
two particles are separated by a (spherical) "monolayer" 
of third particles of the first coordination shell. In order 
to come closer together, the two particles need to diffuse 
"uphill" , i.e. jump over the particle in between. In this 
case the diffusive driving force is pointing outwards along 
r and it generates an additional repulsion with respect to 
the bare local osmotic pressure which is reflected in the 
higher local maximum of V e s as compared to the case 
of hard spheres in vacuo. This effect is particularly im- 
portant for diffusion-limited reaction kinetics in crowded 
systems HHHHI). 
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FIG. 4: (Color online) Schematic representation of the model 
used: an excluded- volume chain in a suspension of Brownian 
particles. The background solvent, taken to be structureless, 
is not shown. 
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III. COIL-GLOBULE TRANSITION DRIVEN 
BY CROWDING 




0.05 0.1 0.15 0.2 0.25 
Volume fraction of crowders <j> 

FIG. 5: (Color online) The values of B2 calculated numer- 
ically from Eq.©, showing the point of changing sign at a 
critical concentration (f> c « 0.145. The straight line illustrates 
the linear-order approximation, B2 ~ (27t<t 3 /3)(1 — (17/4)0o), 
which strongly overestimates the critical point cj>c- 



Let us now imagine that the two interacting particles 
represent two monomers of a polymer chain in a suspen- 
sion of other Brownian molecules (crowders) which in- 
teract with the two test monomers via excluded volume. 
The model is sketched in Fig.|4j For this system, we note 
that the effective interaction derived above is a function 
of the overall density in the system, which includes both 
the monomers and the crowders. In a more precise way 
we could decompose the two contributions such that: 

0o 4>tot 4>tot 

where 4> rnc is the local profile of crowders around the 
test monomer, <p mm is that of other monomers, and 
0tot is the overall bulk packing fraction. At high den- 
sity of crowders, <p mc 3> mm and the general radial 
distribution function is determined by crowders only 
g(r) ?s 4>mc(j , )/(f>o, while in the opposite limit of 4> mc — > 
we have g(r) w 4> m m(r)/(j)o, that is, only determined by 
the hard-sphere interaction of monomers. We can write 
down the effective second virial coefficient B2 for the two 
chosen polymer segments. Normally one expects B2 to be 
a function of T (as well as the mean concentration </>n), 
however, in our case this expression is purely cntropic 
and therefore completely athermal: 



Bo = 2tt 



1 



dr. 



(8) 



The sign of the second virial coefficient controls whether 
the polymer chain would be found in the collapsed state, 
which corresponds to B2 < 0, or swollen (B2 > 0) [29| . 
as we are going to discuss below. For a given </>n the sign 
of B2 is dictated by the balance in the integrand between 
the repulsive hard-core part of the effective potential (for 
spheres of diameter a: B2 = 2-7r<7 3 /3) and the attrac- 
tive part due to the unbalanced osmotic pressure. As the 



hard-core part is independent of 4>q , the behavior of B2 as 
a function of the volume fraction of crowders is therefore 
controlled by the behavior of the attractive part of V c ff(r) 
as a function of </>o- As shown in Fig. [3J the strength of 
the attractive part of V e s (r) is a monotonically increasing 
function of (f>o- This is intuitively expected because the 
unbalanced osmotic pressure in the excluded volume gap 
between two monomers (the standard depiction) and the 
diffusive driving force both become stronger with increas- 
ing density of crowders. Figure [5] shows the evolution of 
B2 with increasing concentration of crowders, together 
with the linear-order approximation valid at very low 
densities, which is easily obtained by the series expansion 
of Eqs.© and ([5]). There is a remarkable similarity with 
the traditional T-dcpcndence of B2(T) for the Lennard- 
Jones type of attractive-repulsive potentials, however, we 
need to emphasize once again that our basic theory is 
purely entropic and therefore athermal. The application 
of the effective potential V e s from Eq.© to evaluate B 2 , 
leads to a change of sign of the second virial coefficient 
in a system of purely hard-core repulsive particles, when 
a critical concentration cf> c 0.145 is reached. 

As a simple exercise to illustrate the theory we take the 
Flory excludcd-volume theory of polymer chains, where 
the free energy is given by the two contributions both 
determined by the size (radius of gyration R) of the 
chain [29| , the total potential energy of effective pair in- 
teractions and the entropy of configurations of an ideal 
(Gaussian) random walk: 



F~kT[B 2 



N 



n 2 



b 2 ') 



(9) 



where N is the number of monomers in the chain and 
b is the linear size of the monomer. To simplify the 
key qualitative results, let us assume b = er, that is, 
the size of the chain segments and the crowders is the 
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same, as illustrated in Fig. [4] In spite of its many con- 
ceptual shortcomings, the Flory expression © is remark- 
ably accurate. Although it neglects correlations between 
monomers along the chain, the errors in the internal en- 
ergy and the entropic terms are roughly the same, which 
leads to the cancelation of errors responsible for the good 
accuracy of this crude model [29[. Minimizing Eq.© 
with respect to R leads to the classical estimate of the 
size of polymer chain swollen in good solvent: 



R ~ {B 2 YN 6 ) 



2xr3U/5 



(10) 



which applies for B 2 > 0. When B 2 < the free en- 
ergy ([9]) does not have a well-defined minimum and the 
chain is considered collapsed into a dense globule. The 
theory of polymer collapse in a poor solvent is developed 
[30H32j . and it is known that to properly describe the 
globular (collapsed) state one needs to invoke the third 
virial coefficient contribution to the interaction energy, 
which provides the stabilizing effect against an apparent 
collapse to a point, R — > 0. In this paper we do not deal 
with the globular state and the considerable computa- 
tional difficulties associated with it. Our aim is merely 
to illustrate the effect of molecular crowding on a chain 
in an otherwise good solvent, and describe the onset of 
the chain collapse, which occurs at B 2 > where B$ is 
not yet relevant. 

Using Eq.© for the effective potential evaluated mak- 
ing use of Percus-Yevick theory for g(r) [26| in Eq.© 
and the latter in Eq.@, the polymer size as a function 
of the crowder concentration c/jq has been evaluated for 
b = a and plotted in Fig. [5] The power-law scaling with 
the crowder density follows directly from Eq. lfTOf when 
one substitutes the linear relation B 2 m const • (<f> c — cf>o) 
near the transition point. 

One may argue that the characteristic scaling of the 
critical point might change with the third virial coeffi- 
cient taken in to account. Therefore, we also made an es- 
timate of the polymer size based on the perturbative ex- 
pansion by Muthukumar and Nickel [34j , which accounts 
for the effect of the three-body collisions and for chain 
connectivity in a simple way. The perturbation is done on 
the distribution function of the end-to-end vector using 
B 2 as the small parameter near the point of its chan ging 
sign. The full expression is reported as Eq.(2.100) in [33( 
where, strictly, the perturbation expansion parameter is: 
z = (3/2n) 3 / 2 B 2 N 1 / 2 /b 3 . A few points calculated at 
(4>o — 4> c ) — > with this approach are shown in Fig. [6] 
and are found to effectively coincide with the more naive 
Flory calculation. Hence we will stick to the handier 
Flory description to analyze the physics predicted by the 
model. 

The theory predicts a sharp collapse of the polymer 
at a well-defined critical value of the crowder volume 
fraction, at cf> c ~ 0.145, which is the point at which B 2 
changes sign and becomes negative. This value is close 
to the result of [23[ where it was experimentally found 
that DNA, under conditions of screened electrostatics, 
collapses sharply when 10-20% of PEG is added to the 
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FIG. 6: (Color online) The size of the polymer as a function 
of the volume fraction of molecular crowders for 6 = a — lnm 
and number of monomers N = 300. The critical transition 
occurs at the point 0n = <t>c where the second virial coefficient 
B2 changes its sign. At high concentration the collapsed poly- 
mer chain is assigned a dense-packed size R ~ WV 1/3 . The 
additional data points (A) have been calculated using the 
perturbation expansion for excluded volume. 



solution. In the limit of infinite dilution of crowders, 
<^o = 0, the effective potential for the monomers reduces 
to the pure hard-core interaction and B 2 = 2nb 3 /3 (since 
we treat the background solvent as a structureless contin- 
uum) and thus we recover the limit of fully swollen self- 
avoiding chain in athermal solvents: R = (2n/3)bN 3 ^ 5 . 
Previous approaches treat the crowders like Newtonian 
hard spheres with a uniform local density and, therefore, 
underestimate the attraction. In this context, the the- 
ory of [2(| predicts no collapse for systems where crow- 
ders and monomers are comparable in size, contrary to 
the evidence of collapse also in that limit from experi- 
ments [H]. 

In the present study we have focused, to preserve an- 
alytical clarity of the key points, on crowding molecules 
which have a size comparable to the size of the monomers 
in the polymer chain. In future studies, this theory 
should be extended to arbitrary monomer/crowder size 
ratios. The theory in fact accounts for this ratio through 
the radial distribution function g(r) of crowders around 
a monomer, which can be extended to the case of arbi- 
trary size ratios by means of a theory by Lebowitz [35| . 
Unfortunately, this can be done only numerically. How- 
ever, we can anticipate that, upon increasing the size 
of the crowders, the packing around the monomer be- 
comes less efficient and the first peak in the g(r) get s 
smaller at the same mean packing fraction (j>o [l|, |35j. 
This means that the collapse transition occurs at higher 
(j)o upon increasing the crowders size with respect to 
the monomer. The reverse occurs upon decreasing the 
crowder/monomer size ratio: the first peak is higher at 
comparatively lower 4>q and therefore the collapse oc- 
curs at lower <fio too. Another interesting direction would 
be to extend this theory to specific practically relevant 
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biomolccular systems by simply adding the specific inter- 
actions (hydrophobic, hydrogen-bonding, etc, as well as 
the chain bending rigidity) that are relevant in each case, 
e.g. in packing the DNA or (mis)folding proteins. These 
additional potentials act as perturbations by shifting the 
critical concentration of the coil-globule collapse, but the 
effect of crowding is still controlled by the effective two- 
particle potential V e fi(r), which is essentially entropic, 
excluded-volume driven. Previous approaches [20j omit 
the microscopic details of the local structure induced by 
the crowders and therefore are not suitable for a realis- 
tic description of the interaction chemistry beyond the 
hard-sphere interaction. 



tcntially capable of accounting for non-trivial chemistry 
details, whereas earlier approaches to depletion cannot 
go beyond the hard-sphere interaction between monomer 
and crowder. The enhanced attraction due to the diffu- 
sive drift in a structured environment of crowders leads 
to the collapse of the polymer, in agreement with simula- 
tion and experimental results, whereas previous studies 
neglecting the diffusive effect (and thus underestimating 
the attraction) predict that collapse cannot occur in that 
limit. 
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IV. CONCLUSIONS 

In summary, we have developed a theory that ana- 
lytically describes the effective interaction between two 
particles in Brownian hard-sphere liquids. In addition to 
the multi-particle correlations which induce the depletion 
attraction in hard-spheres in vacuo (governed by Newto- 
nian dynamics), we also account for the diffusive driving 
force acting on the two Brownian particles due to the 
local gradients of osmotic pressure and chemical poten- 
tial dictated by the liquid-like structure. In fact, when 
two particles move into the depletion gap (near-contact 
region) they are subject not only to the unbalanced os- 
motic pressure in the gap (the cause of the standard de- 
pletion attraction), but also to the additional diffusive 
driving force which pushes them towards low-density re- 
gions (the depiction gap) thus enhancing the depiction 
attraction. Slightly further apart, the particles also ex- 
perience a much higher repulsive barrier where the local 
density decreases between the first and the second coor- 
dination spheres. The effective potential V c g(r) uses the 
structure factor g(r) to calculate the forces between the 
pair of particles mediated by the structured dense fluid 
dominated by Brownian motion. 

It turns out that this additional effect can be very 
substantial in systems of equal-size particles where it in- 
creases the attraction by over a factor of 4 (at 0o = 0.25) 
in comparison to the standard depletion interaction term 
— In g(r) for hard spheres in vacuo. The latter form is cor- 
rectly recovered in the low density limit <j>Q — > 0, where 
the local gradients of density and osmotic pressure dis- 
appear and there is no difference between hard spheres 
in vacuo and Brownian hard spheres. 

As an illustration of the theory, we applied it to cal- 
culate the size of a polymer chain in solution with Brow- 
nian particles (crowders) of size equal or comparable to 
the monomers. Our theory predicts a sharp coil-globule 
transition occurring at a packing fraction of crowders 
(f> c rj 0.145, which is the volume fraction at which the sec- 
ond virial coefficient of the monomers changes sign owing 
to the effective pair potential. This approach is fully mi- 
croscopic because it uses as input the radial distribution 
function of the monomcr/crowdcr system which is po- 
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APPENDIX: ANALYTICAL FORM OF g{r) 

In order to calculate the potential of mean force be- 
tween monomers, we need to evaluate the radial distribu- 
tion function g(r). It is possible to do so analytically for 
the case of spherical crowders which have size comparable 
to the monomers and interact purely by excluded volume. 
In this case the g(r) is the one of hard-sphere fluids at 
given density pa 3 and volume fraction </>o = irpa 3 /6. The 
following analytical representation of g(r) has been used 
for our calculations, which is based on the analytical solu- 
tion of the Percus-Yevick equations by Wertheim [36[ and 
the successive improved paramcterizations by Trokhym- 
chuk, Nczbeda, Jirsak and Henderson [2(| . According to 
the work of j26[, the g(r) can be constructed by split- 
ting it into a depletion part valid for a < r < r* and 
a structural part valid for r > r* . The depletion part 
describes the first coordination shell up to the first (de- 
pletion) minimum whereas the structural part describes 
the damped oscillations at farther distance. The analyt- 
ical expressions read as follows (2(| [37| : 

g(r) = -e^ r -^ + — cos[/3(r - a) + -/]e ai - r ~^ 
r r 

for a < r < r* 

c 

= 1-1 cos[wr + 5}e~ Kr 

r 

for r > r* 

The parameters which figure in the above expressions are 
all ultimately functions of </>o only and can be evaluated 
using the improved parameterizations of (26j based on 
thermodynamic consistency arguments. The parameters 
read as: 

A = <jg a - B cos 7 (12) 

B _ g m -(agjr*)e^'-^ ^ 

cos[/3(r* - a) + ^}e a( - r '-^ - cos[7]e^( r *- CT ) T 

(13) 
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C 



r*{g m -l)e Kr 
cos[o;r* + S] 



(14) 



where 
ft 



6 = -tor* - arctan[(«x* + l)/(ur*)} (15) 



g m =1.0286 - 0.609500 + 3.5781$; - 21. 365105 

+ 42.63440^ - 33.84850^ 
I 



V = 2.0116 - 1.064700 + 0.053800* (18) 



aa = 44.554 + 79.8680 o + 116.4320^ - 44.652e 2<#, ° (19) 



,8a = -5.022 + 5.8570 o + 5.089e~ 4 *° (20) 



7 = arctan 



1 a{a (a 2 + 0g) - ^(a§ - /3g))(l + ±0 O ) + (ag + $ - //«o)(l + 20 o ) 



ft 



a(a 2 Q fxa )(l + ±0o) - n(l + 20 o ) 



(21) 



r 



20o / d 0o 

aoa = T^t 1 + i07"2d 



Furthermore, the parameters k and w have been evalu- 
ated according to the parameterization of Roth et al. [38| 
(22) an d rea< A as follows: 



o 2 o R ( d 



(23) 



ko- = 4.674e- 3 - 93500 + 3.536e~ 56 - 27000 (26) 



20c 



20o + d 



(24) 



= -O.682e~ 24 - 69600 + 4.720 + 4.45O0 O . (27) 



d = [200(0^-300-3+^/3(04 - 20g + 0g + 60o + 3))] 1/3 

(25) 



[1] J. P. Hansen and I.R. McDonald, Theory of Simple Liq- 
uids (Academic Press, London, 2006). 

[2] M. Dijkstra, R. van Roij, R. Evans, Phys. Rev. E 59, 
5744 (1999). 

[3] T. Biben, P. Bladon, and D. Frenkel, J. Phys.: Con- 
densed Matter 8, 10799 (1996); R. Dickman, P. Attard, 
V. Simonian, J. Chem. Phys. 107, 205 (1997); A. Zac- 
cone, H. Wu, and E. Del Gado, Phys. Rev. Lett. 103, 
208301 (2009). 

[4] S. Asakura, F. Oosawa, J. Chem. Phys. 22, 1255 (1954); 
R. Roth, R. Evans, and S. Dietrich, Phys. Rev. E 62, 5360 
(2000). 

[5] D. Marenduzzo, K. Finan, and P. R. Cook, J. Cell. Biol. 
175, 681 (2006). 

[6] J. Dhont, An Introduction to Dynamics of Colloids (El- 
sevier, Amsterdam, 1996). 

[7] L. Lerman, Proc. Natl. Acad. Sci. USA 68, 1886 (1971). 



[8] V. A. Bloomfield, Curr. Opin. Struct. Biol. 6, 334 (1996); 

Biopolymers 44, 269 (1997). 
[9] V. B. Teif and K. Bohinc, Prog. Biophys. Mol. Biol. 105, 

208 (2011). 

[10] L. Tortora and O.D. Lavrentovich Proc. Natl. Acad. Sci. 

USA 108, 5163 (2011). 
[11] T. P. J. Knowles, C. A. Waudby, G. L. Devlin, S. I. A. 

Cohen, A. Aguzzi, M. Vendruscolo, E. M. Terentjev, M. 

E. Welland, and C. M. Dobson Science 326, 1533 (2009). 
[12] M. S. Cheung, D. Klimov, and D. Thirumalai, Proc. Natl. 

Acad. Set. USA 102, 4753 (2005). 
[13] S. Edwards and M. Muthukumar, J. Chem. Phys. 89, 

2435 (1988). 

[14] A. Baumgartner and M. Muthukumar, J. Chem. Phys. 

87, 3082 (1987). 
[15] M.R. Shaw and D. Thirumalai, Phys. Rev. A 44, R4797 

(1991). 



9 



[16] A. Vrij, Pure and Appl. Chem. 48, 471 (1976). 

[17] A. Yu. Grosberg, I. Ya. Erukhimovitch, and E.I. 

Shakhnovitch, Biopolymers 21, 2413 (1982). 
[18] R. P. Sear, Phys. Rev. E 58, 724 (1998). 
[19] P. van der Schoot, Macromolecules 31, 4635 (1998). 
[20] M. Castelnovo and W.M. Gelbart, Macromolecules 37, 

3510 (2004). 

[21] J. S. Kim and I. Szleifer, J. Phys. Chem. C 114, 20864 
(2010). 

[22] Ya. B. Zeldovich, A. A. Ruzmaikin, and D.D. Sokoloff, 
The Almighty Chance (World Scientific, Singapore, 
1990). 

[23] V.V. Vasilevskaya, A. R. Khokhlov, Y. Matsuzawa, and 

K. Yoshikawa, J. Chem. Phys. 102, 6595 (1995). 
[24] N. Dorsaz, C. De Michele, F. Piazza, P. De Los Rios, and 

G. Fofii, Phys. Rev. Lett. 105, 120601 (2010). 
[25] J.L. Barrat, J. P. Hansen Basic concepts for simple and 

complex liquids (Cambridge Univ. Press, Cambridge, 

2003). 

[26] A. Trokhymchuk et at, J. Chem. Phys. 123, 024501 
(2005). 

[27] J. Dzubiella, J.A. McCammon J. Chem. Phys. 122, 



184902 (2005). 

[28] A. Zaccone, N. Dorsaz, F. Piazza, C. De Michele, M. 
Morbidelli, and G. Foffi, J. Phys. Chem. B 115, 7383 
(2011). 

[29] M. Rubinstein and R. H. Colby, Polymer Physics (Oxford 

University Press, Oxford, 2003). 
[30] TA. Vilgis Phys. Rep. 336, 1677254 (2000). 
[31] A.L. Kholodenko and K.F. Freed KF J. Phys. A: Math. 

Gen. 17 2703 (1984). 
[32] V. G. Rostiashvili, G. Migliorini, and TA. Vilgis, Phys. 

Rev. E 64, 051112 (2001). 
[33] M. Doi, S.F. Edwards The Theory of Polymer Dynamics 

(Clarendon Press, Oxford, 1986) . 
[34] M. Muthukumar and B.G. Nickel, J. Chem. Phys. 80, 

5839-5850 (1984). 
[35] J.L. Lebowitz Phys. Rev. 133, 895-899 (1964). 
[36] M.S. Wertheim, Phys. Rev. Lett. 10, 321 (1963). 
[37] A. Trokhymchuk et al., J. Chem. Phys. 124, 149902 

(2006). 

[38] R. Roth, R. Evans, and S. Dietrich, Phys. Rev. E 62, 
5360 (2000). 



